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ABSTRACT 

It has long been recognised that quasi-periodic oscillations (QPOs) in the X-ray light curves of accreting black 
hole and neutron star binaries have the potential to be powerful diagnostics of strong field gravity. However, 
this potential cannot be fulfilled without a working theoretical model, which has remained elusive. Perhaps the 
most promising model associates the QPO with Lense-Thirring precession of the inner accretion flow, with the 
changes in viewing angle and Doppler boosting modulating the flux over the course of a precession cycle. Here, 
we consider the polarization signature of a precessing inner accretion flow. We use simple assumptions about 
the Comptonization process generating the emitted spectrum and take all relativistic effects into account, paral¬ 
lel transporting polarization vectors towards the observer along null geodesics in the Kerr metric. We find that 
both the degree of linear polarization and the polarization angle should be modulated on the QPO frequency. 

We calculate the predicted absolute rms variability amplitude of the polarization degree and angle for a specific 
model geometry. We find that it should be possible to detect these modulations for a reasonable fraction of 
parameter space with a future X-ray polarimeter such as NASA’s Polarization Spectroscopic Telescope Array 
(PolSTAR: the satellite incarnation of the recent balloon experiment X-Calibur). 

Subject headings: polarization - accretion, accretion disks - black hole physics - X-rays: binaries 


1. INTRODUCTION 

Low frequency quasi-periodic oscillations (hereafter 
QPOs) with frequencies ranging from ~ 0.1 — 30 Hz are 
routinely observed in the X-ray light curves of black hole 
(BH) and neutron star binary systems. The QPO proper¬ 
ties correlate strongly with spectral state, which evolves 
from the hard power law dominated hard state to the disk 
blackbody dominated s oft st ate on timescales o f ^months 
(iTananbaum et al.l 1 19721 see Ivan der Klisl 120061: iDone et all 
120071 iBellonil 120101 for reviews). The QPO frequency 
increases as the spectrum softens and the amplitude r ises t o 
a peak in the hard intermediate state (HIMS: IBellonil 1201 Oh 
before reducing as spectral evolution f urther continues (e.g. 
iMuno et ahl 1 1 9991 ISobczak et al.ll20()0h . Since the flux also 
peaks in this state, the HIMS is ideal for studying QPOs. 

The disk blackbody spectral component is from a ge¬ 
ometrically thin accretio n disk (IShakura & Sunvaevl 119731 
iNovikov & Thornall973l) and the power law component re¬ 
sults from Compton up-scattering of cool seed photon s by 
some cloud of hot electrons close to the BH ( iThorne & Price! 
119751 ISunvaev & Truemperll 1 97 c jll . Some fraction of these 
seed photons are provided by the disk, with the rest gen¬ 
erated inte rnally in the electron cl oud via cyclo-synchotron 
radiation (iGhisellini et alJ 119881 iPoutanen & Vurml 120091 
IVeledina et al.l201 lbl) . The exact geometry of the Comptoniz- 
ing cloud is still a matter of debat e but a prominent interpre¬ 


tation is the truncated disk model (Ichimaru 

,1977||Esin et alJ 

IT997^|Poutanen et alTl997: Done et al. 2007 

; Gilfanovf2010). 


in which the geometrically thin disk evaporates inside of some 
radius larger than the innermost stable circular orbit (ISCO) 
to form a geometrically thick, optically thin (r ~ 1) accre¬ 
tion flow (hereafter the inner hot flow). The spectral transi¬ 
tions can be explained if the truncation radius moves smoothly 
from R 0 ~ 60 R g , where R g = GM/c 2 , in the hard state to 


R 0 = -Risco in the soft state. 

It has long been recognised that QPOs have the poten¬ 
tial to be powerful diagnostics of the regions close to BHs. 
However, this potential cannot be realised without a work¬ 
ing model, which has long remained elusive. Suggested QPO 
mechanisms in the litera ture either consider particle orbits in 
Gener al Relativity (GR) (I Stella & Vietrilfl998l : IWagoner et alJ 
120011: ISchnittman et alJ l2006al) o r instabilities in the accre¬ 
tion flow dTagger & Pellatl 119991 : ICabanac et alJl20Tol) . Per¬ 
haps the most successful model for explaining the array of 
observational properties attributes the QPO to the effect of 
frame dragging. In GR, a spinning massive object twists up 
the surrounding space-time, leading to precession of particle 
orbits inclined with the BH equatorial plane. This is gener¬ 
ally called Lense-Th irring precession after t he authors who 
originally derived it (iLense & Thirringll 19181) , although their 
derivation was only in a weak fie ld limit, coming well before 
the de rivation of the Kerr metric (iKerril 19631) . IStella & Vietril 
(119981) first suggested that the QPO results from Lense- 
Thirring precession, considering only the p recess ion fre¬ 
quency of test masses at different radii. ISchnittman et alJ 
( 2006a;) considered instead a precessing ring in the accre¬ 
tion disk. However, phase resolved spectroscopy reveals that 
it is the Comptoni zed spectral componen t which oscillates 
and n o t the dis k (Markwardt - etal J 1 1 9991: iR evnivtse vetal J 
20011: ISobolewska & Zvckil 120061 : lAxelsson et al.l 120131) . 
Ingram et alJ (120091) suggested that the entire inner flow 
precesses as a solid body, motivated by the General Rel¬ 
ativistic Mag neto-Hydrodynamic (GRMHD) simulations of 
[Fragile et al.l ( 2007 ). This precessing flow model naturally 
explains the observed QPO s pectrum, pre dicts the correct 
range of frequencies for BHs (llngram et al.ll2009l) . has been 
incorporated i nto a full model for the power spectral prop- 
erties of BHs ( fingram & Done 2011 : Ingram & Donejf2012al : 
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Ilngram & van der Klisl 1201 3t iRanisarda et al.ll20l4l) and has 
been exten ded to also account for the optical QPOs observe d 
from BHs (IVeledina et alJl20l3t IVeledina & Poutanenl[2~015l) . 
The model can also explain the range of QPO frequencies 
observed in low accretion rate neutron stars (Ilngram & Donel 
1201 Oh. but not the 11 Hz p ulsar in the globular cluster Terzan 
5 (lAltamirano et aT]l2012f) . 

In this Paper, we consider the X-ray polarization sig¬ 
nature from a processing inner flow. Photons emerging 
from a Comptonizing slab are expected to be polarized, 
with the polarization degree as a function of viewing an¬ 
gle dependingjm o ptica l depth of the sla b dCha ndrasekhari 
1960tlRvbicki & LightmaJl 9791: [Loskutov & Sobolev! 19821 : 

Sunvaev & Titarchukf 1985ciPoutanen & Svenssonll 19%[h Al¬ 

though the polarization degree is predicted to be less than 
~ 20%, the current generation of proposed X-ray polarime- 
try missions should comfortably be able to detect a sig¬ 
nal. Polarimetry promises to provide a powerful extra 


lever arm to interpret the 
predicted energy dependen 
has been extensively explot 

X-ray signal from BHs. The 
ce of the polarization signature 

ed (e.g. Stark & Connors 1977; 

Dovciak et all 2008b! iLi e 

all 2009; Schnittman & Krolik 

2009; Schnittman & Krolik 

201(| Krawczvnskll2012). Vari- 


ability of the signal is less well studied, although Lense- 
Thirring precession did provide an early motivation for X- 
ray polarimetery in the 1970s (Kn ox Lo ng, p rivate com- 
munication). More rec ently, iDovciak et al.l (l2008al) and 
iZamaninasab et all (120111) co nsidered the polarization signa- 
ture of orbiting hotspots and iMarin & Dovciakl (120151) con¬ 
sidered the polarization signatures of obscuring clouds, all for 
the case of Active Galactic Nuclei. In this Paper, we predict 
a modulation in the polarization signature from bright stellar 
mass BHs on a timescale of ~ Is. We use simple parame- 
terizations for the angular dependence of emissi on and polar¬ 
ization designed to mimic the results of lSunvaev & Titarchukl 
(fl985l hereafter ST85) for the r = 1 case. We then take 
GR fully into account by ray tracing photon paths from the 
precessing flow geometry to the observer, assuming the Kerr 
metric. 

In Sections [2] and [3] we outline our assumed geometry and 
formalism for calculating the observed polarization degree 
and angle as a function of QPO phase. In Section^ we con¬ 
sider a specific flow geometry appropriate for the HIMS and 
a QPO frequency of ~ 1 — 2 Hz. We calculate the ampli¬ 
tude of the modulation in both polarization degree and angle 
for the full range of observer positions. In Section [5] we dis¬ 
cuss the prospects of detecting this signature and suggest im¬ 
provements that can be made to our modelling assumptions 
in future. Throughout this Paper, we express distance in units 
of R g and time in units of c/R g in order to exploit the scale 
invariant nature of GR. We use Einstein’s summing conven¬ 
tion with Greek indices taken to run from 0 to 3 and Roman 
indices taken to run from 1 to 3. We adopt the 4-vector for¬ 
malism with negative time entries and positive spatial entries. 

2. GEOMETRY 

We ass ume the geometry propos ed b vllngram et alj(l2009l). 
following Ilngram & Donel (12012bl) and IVeledina et al.l (120131 
hereafter VPI13). Specifically, the BH and binary system spin 
axes are misaligned by some modest angle /3. The outer thin 
disk is assumed to be in the binary plane and the flow spin axis 
is assumed to precess around the BH spin axis, always main¬ 
taining the same misalignment angle /?. In this section, we 
first describe our coordinate system followed by a description 



Figure 1. Coordinate system used in this Paper. We define right handed 
Cartesian coordinate systems with z-axes aligned with the BH. flow and bi¬ 
nary spin axes. The BH (black) and binary (red) z-axes are misaligned by 
the angle 0 and the binary system x-axis is chosen to be in the unique plane 
shared by the BH and binary z-axes. The vector 6 represents the observer’s 
line of sight, and the BH x-axis is chosen as the projection of this vector onto 
the BH equatorial plane. The flow z-axis (blue) precesses around the BH 
z-axis with precession angle tu, always maintaining the misalignment angle 
0. The flow and binary z-axes are therefore aligned when ui = 180° and 
misaligned by an angle 20 when ui = 0°. 

of the assumed geometry of the inner flow. 

2.1. Coordinate system 

FigurejTjillustrates the coordinate system used in this paper. 
As in VPI13, we assume the binary spin axis is misaligned 
with the BH spin axis by an angle /3. We define a ‘binary’ 
coordinate system in which the z-axis, Zb, aligns with the bi¬ 
nary spin axis and therefore the plane of the binary is simply 
the plane perpendicular to Zb- We choose to align the binary 
system x-axis with the projection of the BH spin axis on the 
plane of the binary. In binary coordinates, the vector pointing 
from the BH to the distant observer is given by 

6 = (sini cos$, sin* sin4>, cos*), (1) 

where i is the inclination angle and <1> is the viewer azimuth. 
We use these angles to define the observer’s position be- 
cause it is possib l e to m easure i via dynamical methods (e.g. 
lOrosz et al J2004l201 ll) . Throughout this Paper, a hat denotes 
a unit vector. 

We assume that the flow spin axis (Zf) precesses around 
the BH spin axis, maintaining a misalignment angle /3 as the 
precession angle uj increases. This means that the misalign¬ 
ment between the flow and binary spin axes varies over a pre¬ 
cession cycle between 0 and 2/3. We define a ‘flow’ coordi¬ 
nate system with exactly the same relationship to the binary 
coordinate system as introduced in VPI13. Here, however, 
we will perform calculations in the Kerr metric which is az- 
imuthally symmetric around the BH spin axis, in contrast to 
the Schwarzschild metric used in VPI13 which is spherically 
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symmetric. Consequently, our calculations are greatly simpli¬ 
fied by defining a ‘BH’ coordinate system with basis vectors 
x, y and z. The z-axis aligns with the BH spin axis, allowing 
us to use Boyer-Lindquist coordinates, and the x-axis aligns 
with the projection of 6 onto the BH equatorial plane. In this 
coordinate system, the observer’s line of sight can be written 
as 

6 = (sin0 o ,O,cos0o)i (2) 

where 

cos0q = sin i cos <f> sin/3 + cos i cos/3, (3) 

is the cosine of angle between the observer’s line of sight and 
the BH spin axis. We can express the flow basis vectors in this 
coordinate system as 

Xf = (— cos/3cos(w — wo), — cos/3sin(w — wo), sin/3) 
y f = (sin(w - w 0 ), - cos(w - w 0 ), 0) 

Zf = (sin/3cos(w — wo), sin/?sin(w — wo), cos /3), (4) 

where w>o is the precession angle at which the projection of Zf 
onto the BH equatorial plane aligns with the BH x-axis (see 
Figure©. This is given by 

sin i sin $ 

tanw 0 = ^- . . „ ■ (5) 

sin % cos <f> cos p — cos i sin p 

Note that the flow y-axis always remains in the BH equato¬ 
rial plane and the flow x-axis maintains a constant misalign¬ 
ment with the BH equatorial plane of /3. When w = 0°, the 
flow is maximally misaligned (2/3) with the disk, and when 
w = 180°, the flow is aligned with the disk. Any point on the 
surface of the flow can then be represented in this coordinate 
system as 

r = r (sm0fcos4>f : sm9fsin(/)f,cos9f), (6) 

where 9f and Of are respectively the polar and azimuthal an¬ 
gles defined in the flow coordinate system. The same point 
can alternatively be represented in terms of the BH polar and 
azimuthal angles 

r = r (sin(?cos<)>, sin$sin</), cos$). (7) 

We present the equations to convert between these two sets of 
coordinates in Appendix |C] We can also represent a point in 
terms of binary polar and azimuthal angles, 9b and fa. 

2.2. Inner hot flow geometry 

We assume that the inner flow is shaped like a torus de¬ 
scribed in Boyer-Lindquist coordinates by the equation: 

T(r,9 f ) = B - sin Of [1 - (1 - B)r/r a \, (8) 

where B = [(h/r ) 2 + l] -1 / 2 , r 0 is the outer edge and h/r is 
the scale height of the flow. Outside of the flow, T(r,9f) > 0 
and inside the flow, T(r,9f) < 0. On the flow boundary, 
T(r,9f) = 0. Figure [2] shows some example cross-sections 
of this shape; i.e. contours of T = 0. For r -C r D , the flow 
has a constant opening angle, with sin 9f = B. For larger r, 
sin 9 f > B and the function curves round in the shape of a 
tear drop. We also assume an inner radiation edge, r;, which 
we set equal to the ISCO in this Paper. We can use Equation 
([8]) to assess if a given position, described by the vector r, is 
inside, outside or on the boundary of the flow. 

This shape is motivated in part by the discussion on thick 
disks in Chapter 10 of iFrank et akl (120021) . Equation © 
comes from considering an inviscid fluid rotating about a 



rameters: the scale height h/r and the outer radius in units of R g , r 0 . For 
r -C r 0 , this function has a constant opening angle, but curves around into a 
tear drop shape for larger r. We show three example parameter combinations: 
h/r = 0.1, r„ = 20 (solid line), h/r = 0.2, r 0 = 25 (dashed line) and 
h/r = 0.4, r 0 = 30 (dot-dashed line). In this Paper, we consider the first of 
these geometries, represented by the solid line. 

central (Newtonian) point mass (see Figure 10.2 therein). 
Clearly this is over simplified, but it does provide a conve¬ 
nient way to parameterize the flow geometry in the most re¬ 
alistic way currently possible. Indeed, the shape illustrated 
in Figure [2 is comparable to that seen in GRMHD simula¬ 
tions (e.g. lFragile et alj|2007t IFragile & Meietll2009l: iFragilel 
120091) which were initialized with a ‘Polish doughnut’ so¬ 
lution (iJaroszvnski et all 1 1 9801) and allow ed to evolve self- 
consis tently to show solid body precession. iDexter & Fragilel 
(1201 ll) used ‘after the fact’ assu mptions about the radiat ive 
emissivity of the flow (following lSchnittm an et al l12006b!) in 
order to ray trace emission from the IFragile et al .1 (120071) sim- 
ulation but were unable to study long enough time scales to 
see the precession period modulate the light curve (they con¬ 
centrated instead on trying to find high frequency QPO candi¬ 
dates). Thus defining a reasonable flow geometry analytically 
is currently the best way to study the effects of precession on 
the observed emission. 

3. FORMALISM 

In this Section, we describe our method for calculating the 
flux and polarization properties observed from a precessing 
flow as a function of precession angle. Since we only con¬ 
sider linear polarization, the polarization of the signal can be 
described entirely by the polarization degree, p, and angle, \- 

3 . 1 . Inner hot flow properties 

The specific intensity emitted from a patch of the flow sur¬ 
face is, in general, a function of position, photon energy and 
emission angle, 9 e . Our calculations are greatly simplified by 
assuming that these variables can be separated, such that 

I Ee (E e ,r, p e ) = I Eb q e {r) I(p e ), (9) 

where I Ee is the emitted spectrum, q e (r) and l(p e ) are re¬ 
spectively the radial and angular emissivity profiles and p e = 
cos 9 e . In the absence of a standard radial emissivity law 
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Figure 3. Intensity as a function of emission angle assumed for this Paper. 
We use the analytic function I(fi e ) oc 1 + 1.1(1 — fi e ) — 1.4(1 — /x e ) 2 , 
designed to mimic the results of ST85 for the r = 1 case. This Figure can be 
compared with Figure 4 of ST85, except we use a different normalization. 



Figure 4. Polarization degree as a function of emission angle assumed 
for this Paper. We use the analytic function p(/i e ) = 50%(1 — fi e ) — 
23%[exp {(1 — Me) 1 ' 6 } — 1], designed to mimic the results of ST85 for the 
ro = 1 case. This Figure can be compared with Figure 5 of ST85. 


for a large scale height accretio n flow, we simply assume 
the standard case for a th in disk dShakura & Sunvaevlll973t 
iNovikov & Thorndl 19731) 

q e {r) = r -3 (l - y/ri/r), (10) 

following VPI13. As for the angular profile, we define an 
analytic function designed to mimic the shape obtained from 
the calculations of ST85. We plot this profile in Figure [3 to 
be compared with the tq = 1 case from Figure 4 of ST85. 


Here, we have normalised the intensity to give 

27 T I(fl e )fl e dfl e = 1. (11) 

Jo 

We also parameterize polarization degree as a function 
of emission angle pin,) (defined as the fraction of photons 
emerging from the flow which are polarized) following the 
calculations of ST85. Figure Q] shows this function, to be 
compared with the tq = 1 case in Figure 5 of ST85. We plot 
—p(/i e ) in order to account for the different sign convention in 
ST85. We note that the emission angle 9 e is defined here with 
respect to the flow spin axis rather than to the local normal 
to the torus surface. Thi s is because the calculations of ST85 
(also see IViironen & Poutanenll2004l) assume Thomson scat¬ 
tering from a thin slab and the functions in Figures © and 
0 are therefore only defined with respect to the flow spin 
axis. Since we expect the flow to have a rather small scale 
height ( h/r ~ 0.1) this is likely a fairly good assumption. 
There is however scope to imp rove upon thes e assumpti ons in 
future work. For instance, jPoutanen & Svenssonl (119961) con¬ 
sidered cases where exact solutions can be found outside of 
the Thomson scattering limit. 

3.2. Ray tracing 

We use the publicly available code GEOKERR, described in 

iDexter & Agoll (120091) . to solve for photon geodesics in the 

Kerr metric. We start off by defining an observer’s cam¬ 

era some large distance, D, from the BI-Q along the vec¬ 

tor 6. The impact parameters at infinity, o: ( j and 3q, rep¬ 
resent respectively horizontal and vertical distance on the 

plane of the observer’s camera (a schematic illustrating the 
impact parame ters can be found to the left of Figure Al in 
iMiddleton & Ingraml2015l) . For a given observer position and 
BH spin, a, null geodesics can be uniquely defined by these 
two impact parameters. Alternatively, these geodesics can be 
parameterised by Carter’s constants of motion l = —oto sin 9q 

and q * 2 = (3 q + cos 2 6 *o(q;q — a 2 ). Each combination of im¬ 
pact parameters represents a pixel on the observer’s camera, 
which is hit by a photon on a unique geodesic path. We define 
a grid of impac t parameters with equal logarithmic steps in 

b = y/a.Q + and equal linear steps in the angle ip, defined 
as t&n ip = ao/Po- Ignoring parallax, which is the same for 
each pixel, the solid angle subtended by each pixel is b dbd<p. 

For each pixel we set GEOKERR to compute 100 steps along 
the photon geodesic towards the BH. For each step, a posi¬ 
tion 4-vector x 11 = (t. r, 9, <j>) is provided in Boyer-Lindquist 
coordinates. We convert the BH polar and azimuthal angles 

9 and <j> to the corr espo nding flow angles Of and <l>f using 
Equations (IC3b and ( lC4b . This allows us to evaluate the func¬ 
tion T(r,9f) from Equation ®. If T passes from positive 
to negative between consecutive steps, we conclude that the 
geodesic has crossed the flow boundary. If this does not hap¬ 
pen, we also check if Hf — cos 9f has passed from positive 
to negative, in which case the geodesic has crossed the flow 
mid-plane and therefore must have crossed into and out of the 
flow in a single step. This is most likely to happen for small r 
where the flow is very shallow and a step does not need to be 
particularly large to pass completely through the flow. 

In either of these cases, we have isolated a root of Equation 
(0 between two points on the geodesic path. We use linear 
interpolation to represent rasa function of /// between these 

1 The camera must be far enough away for all geodesics to be straight and 
parallel. We use D = 10 5 * * * 9 * R g . 
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two points, allowing us to express T as a function of Hf only. 
We then find the root of the equation T(fif) = 0 with a bi¬ 
section search. From this solution for fif, we can interpolate 
a solution for r and also <j)f. If r < n, we carry on follow¬ 
ing the geodesic in case it loops back around to intercept the 
underside of the flow. Otherwise, we calculate the contribu¬ 
tion to the flux observed from that pixel. We ignore direct 
disk emission and also emission from the flow reflected from 
the outer disk. This assumption is appropriate for the energy 
range ~ 10 — 20 keV, which is dominated by Comptonized 
emission from the flow. Direct disk emission contributes only 
in soft X-rays < 5 keV and reflection is important at the iron 
line (~ 6.4 keV) and above ~ 20 keV. 

3.3. Disk shielding 

When tracing rays backwards from the observer, we also 
test whether they intercept the outer disk before hitting the 
flow. In this case, since the disk is optically thick, we con¬ 
clude that our view of the flow is blocked for this pixel. We 
assume that the outer disk occupies the binary plane and trun¬ 
cates at r = r Q . At every step along the ray, we convert to 
binary coordinates to assess if fib = cos 9b passes from posi¬ 
tive to negative, indicating a crossing of the disk plane. If this 
happens, we use linear interpolation to calculate the value of 
r at fib = 0. If r > r 0 , we conclude the ray has crossed the 
outer disk and stop following the ray. Otherwise, we carry 
on following it. In the cases where the ray crosses both the 
disk and flow in one step, we use linear interpolation to assess 
which it hit first. 


3.4. Blueshift and flux 

The energy of a photon reaching the observer will be mod¬ 
ified by the gravitational field of the BH and the motion of 
the emitting particle. For a stationary observer at infinity, the 
ratio of the observed to emitted energy of a photon (hereafter 
the blueshift) is given by 


= Eo = -Po 
S E e PeUfj, : 


( 12 ) 


where pi'; and u , ' are respectively the 4-momentum of the pho¬ 
ton and the 4-velocity of the emitting particle (e.g. Luminetl 
119791) . We use a form for the 4-momentum (given in full 
in the Appendix) normalised such that p* = 1, simplifying 
Equation (ITU further. Since we use the metric convention 
of negative time entries and positive distance entries, the 4- 
velocity of the emitting particle is ( u M ) 2 = — (dx^/ds) 2 , 
where ds 2 = g „„ dx M dx v is the line element and g is the 
metric. We can represent this in terms of coordinate veloci¬ 
ties, = dx^ Idt, in the form 

or 

u a = / (13) 

We assume circular orbits (Q r = 0) with Keplerian angular 
velocity about the flow spin axis, so 

d(j)f 1 


n<h = —-I- =_ 

dt r 3 / 2 


(14) 


and fl e f = 0. Note, we assume Keplerian angular velocity 
perpendicular to the flow spin axis even for parts of the flow 
not in the mid-plane. Since VE = 2nv qpo , we convert 


this to BH coordinates as follows: 


dcjif 

de t 
d(j) f L 




(15) 


whe re th e differentials can be computed from Equations (IC5b 
and ( 1C6I >. We use the Kerr metric throughout, which we pro¬ 
vide in Appendix lAl for completeness. 

The specific flux observed at ener gy E „ from a pixel with 
solid angle b d b dp is given by (iLuminedl 19791) 


dF Eo = g 3 lE e (E e , r, fi e )b db dtp. (16) 


For energies significantly greater than the seed photon tem¬ 
perature (fcTbb) and less than the electron temperature (fcT e ), 
the emitted spectrum can be well approximated by a power 
law, I Ee oc El . In this case, the observed spectrum is also 
a power law, Ie oc E^~ Y . Substituting this (and Equation^ 
into Equation (1 161 ). we can represent the flux in a power law 
dominated energy band (~ 10 — 20 keV) as 

dF = g 2+r I Eo q e (r)I(fi e )bdbdUp. (17) 

The emission angle is given by (e.g. iMisner et al.l Il97~3 ': 

iRvbicki & Bromlevlll997HDovciakf20()4l) 


Pe = gPe n pn (l 8 ) 

where n 1 ' is a 4-vector normal to the flow mid-plane, defined 
in the flow rest-frame. See Appendix [B] for the full form of 
this 4-vector. The above formula encapsulates the effects of 
both light bending and relativistic aberration. From this, we 
evaluate dF for every pixel and sum to obtain the total ob¬ 
served flux, F, as a function of precession angle, w. 


3.5. Polarization 

For each geodesic that intercepts the flow, we calculate the 
polarization degree from the function shown in Figure [4] We 
initialize the polarization 4-vector, f 1 ', as the projection of 
n^ on the plane perpendicular Jo the eme rgent photon’s 4- 
momentum, p 1 '. Following. iDovciak ( 20041) . this is defined 
as 


r = 


^ - flejgP^ ~ 

V 1 - El 


(19) 


For each geodesic, we parallel transport the polarization 4- 
vector forwards from the emission point to the observer’s 
camera in geodesic steps, Sx That is, moving from s to 
s + Ss changes the polarization 4-vector to 


r(s+ss) = r(s)-r(s)T : 

with the Christoffel symbols given by 




r CT = -q° 

HV r\ J 


dg VK . , dg, 




dg, 




dxv dx v dx K 


( 20 ) 


( 21 ) 


The position 4-vector, x^, is specified for each point on each 
geodesic by GEOKERR and so we can simply calculate the 
step length as dx 11 = x^(s + Ss) — aF(s). 

Once we have transported f 11 all the way back to the ob¬ 
server, we convert back to BH Cartesian coordinates (see Ap¬ 
pendix 0 and find the projection of f 1 ' on the horizontal and 
vertical directions of the camera. These directions are given 
in BH (Cartesian) coordinates by 

a 0 = (0,1,0) 

do = (— cos 00 , 0 , sin 0 O ), 


(22) 
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Figure 5. Images of the flow for three precession angles for the model with i = 70°, <1? = 110°, /3 = 10° and h/r = 0.1. The left hand 
side pictures blueshift and the right hand side pictures flux with the polarization vector overlaid, normalised to the maximum observed polarization de¬ 
gree. The three precession angles pictured, in units of cycles, are lj = 0, 0.3125 and 0.625 from top to bottom respectively. The full movies for 
these images can be found at http://figshare.com/articles/Polarization_modulation_gifs/1351920 . Indivdual movies can alternatively be found on YouTube: 
www.youtube.com/watch?v=Q2CwOGKVC9U&feature=youtu.be (blueshift) and www.youtube.com/watch?v=E3kYAnS3pQI&feature=youtu.be (flux and po¬ 
larization). 
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and the projections are simply f • d 0 and f • Po, where f is the 
3-vector describing the spatial part of /% The polarization 
angle for a given geodesic is then given by 

f A 

tan[x(a 0 ,/3o)] = —j—A (23) 

I • Po 

We define Stokes parameters for each geodesic 

dQ = dF(a o ,/3 0 ) p{a o ,/3 0 ) cos[2x(ct 0 , A))] (24) 

dU = dF(a o ,/3 0 ) p(a o ,/3 0 ) sin[2x(a 0 , /?o)], (25) 

and sum over all geodesic paths to get the observed Stokes 
parameters Q = f dQ and U = f dU. Note that, under this 
convention, vertical and horizontal polarization corresponds 
to Q/F = 1 and Q/F = —1 respectively for a 100% polar¬ 
ized signal. The overall polarization degree is then 


P = 


\JQ 2 + U 2 


(26) 


and the overall polarization angle, x, can be calculated from 
tan(2x) = k. (27) 

3.6. Simple Newtonian approximation 


In addition to the fully relativistic treatment described 
above, we also define a simple Newtonian approximation to 
assess the importance of GR effects. We make the simplifying 
assumption that the flow is a flat disk, and so any point on the 
flow surface is described by the vector r = xxf + yyt ■ a 
flat metric, photons emitted from r hit the observer’s camera 
at the coordinates ao = r • So and po = r • 0q. The Kep- 
lerian velocity (in units of c) has magnitu de v = r -1 / 2 and 
direction v = (—j/Xf + xyr)/\Jx 1 + y 2 . We use Equation 
(ITTb to calculate the flux for each pixel of the camera. Using 
the Minkowski metric instead of the Kerr metric for Equations 
( IT3] i and i fT2l i. the blue shift becomes 


Vl — v 2 

1 — 6 ■ v 


(28) 


This is the special relativistic Doppler factor. Since all the 
photons emitted at a given time now have the same polariza¬ 
tion angle, the observed polarization vector f can be calcu¬ 
lated as the projection of Zf onto the plane perpendicular to 6 
(if relativistic aberration is ignored). This gives for the polar¬ 
ization angle 


tanx = -- 


sin j3 sin(w — cuq) 


(29) 


sin d 0 cos P — cos f? 0 sin p cos(w — iv njj_ 

which is the standard formula for a rotating vector (IFergusonl 
119731 : IFerguson l: 1976 . IViironen & Poutaneiil 120041) . More¬ 
over, the observed polarization degree simply becomes a flux 
weighted average of the contribution from each pixel. 


4. RESULTS 

In our model, the geometry of the flow is governed by three 
parameters: the inner radius r„ the outer radius r a and the 
scale height h/r. In this Paper, we consider parameters ap¬ 
prop riate for the HIMS and so fix r n = 20 and h/ r = 0.1 
(e.g. Ilngram & Donell201 It llngram & Don el [2012al) . We set 
r'i = nsco and assume a = 0.98. This assumption of high 
spin maximizes the impact of relativistic effects which tend to 
wash out variability, therefore yielding conservative estimates 
for the QPO amplitudes predicted by the model. This corre¬ 
sponds to a QPO frequency of ~ 2 Hz, typical of the HIMS. 


i=70°, 0=110°, (3=10° 



0 0.5 1 


QPO cycles 

Figure 6. Flux, polarization degree and polarization angle plotted against 
QPO phase for the high inclination model. The solid lines are for the fully 
relativistic model and dashed lines representing a simple Newtonian approx¬ 
imation are included for comparison. The fractional rms (for the full model) 
in the first and second harmonics respectively is 10.0% and 2.24% for the 
flux. The absolute rms in the first and second harmonics is respectively 1.0% 
and 0.3% for the polarization degree and respectively 3.1° and 0.2° for the 
polarization angle. 

We additionally fix the misalignment angle to P = 10° fol¬ 
lowing VPI13 and assume a spectral index of T = 2, again 
typical of the HIMS. 

4.1. High and low inclination examples 

We consider the observational appearance of the specific 
geometrical setup described above for a range of viewing an¬ 
gles. Since our geometry is asymmetric, the viewer must be 
specified by both a polar angle i and an azimuthal angle 4>. In 
this section, we first explore two specific examples, chosen to 
represent respectively high and low inclination sources. For 
the high inclination model, we use i = 70° and 4) = 110°. 
For the low inclination model, we use i = 30° and 4? = 180°. 
We then move on to explore a grid of viewing angles in the 
following subsection. For these specific examples, we use a 
high resolution with 400 steps in the impact parameters b and 
ip (i.e. 400 x 400 pixels) and we consider 32 precession an¬ 
gles. 

Figure [5] shows images of the flow for three values of the 
precession angle. The right hand images depict flux, with 
colors defined by the key beneath in arbitrary units. We see 
the characteristic apparent warping of the flow through light 
bending, with the back of the flow appearing to bend towards 
us. Since the BH is spinning rapidly, we also see an asymme- 
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i=30°, <5=180°, (3=10° 



_i_i_i_i_I_i_i_i_i_li 

0 0.5 1 


QPO cycles 

Figure 7. Flux, polarization degree and polarization angle plotted against 
QPO phase for the low inclination model. Solid and dashed lines again re¬ 
fer to the fully relativistic model and the Newtonian comparison respectively. 
The fractional rms (for the full model) in the first and second harmonics re¬ 
spectively is 4.6% and 0.25% for the flux. The absolute rms in the first and 
second harmonics is respectively 1.1% and 0.05% for the polarization degree 
and respectively 10.0° and 1.0° for the polarization angle. 

try to this apparent warping from the frame dragging effect. 
In the top and bottom plots, we see the underside of the flow 
due to rays bending dramatically from their starting point be¬ 
fore passing through the gap between the flow and the disk 
and towards the observer. In reality, we may expect material 
in some transition region between the two accretion regimes 
to block our view of this secondary image to some extent. In 
these images, precession and rotation are anti-clockwise. We 
see that the brightest patch of the flow is always just to the left 
of the BH. This region not only has the largest radial emissiv- 
ity, but emission from here is Doppler boosted by the rapid 
rotation of material in the flow. Also, parts of the flow are 
blocked from view by the outer disk. This can be seen for all 
three precession angles but is perhaps clearest in the top plot 
where the bottom left corner of the flow is hidden behind a 
disk which is not pictured. 

The images on the left show blueshift for the same three 
precession angles, with a key again included beneath. Emis¬ 
sion from the approaching material to the left of the BH is 
blue shifted by Doppler effects, whereas the receding mate¬ 
rial to the right of the BH is red shifted. Close to the BH, 
we see the effects of gravitational red shift, whereby photons 
lose a significant amount on energy escaping the gravitational 
pull of the BH. For selected pixels on the right hand plot, we 
also represent the polarization vector with a black line. The 


length of the line gives the magnitude of the polarization de¬ 
gree, normalised to the maximum measured from the entire 
run. The orientation of these v ectors is heavily influen ced by 
GR through two main effects (IStark & Connorslil97? ). First 
of all, light bending means that photons reaching the observer 
may have had a different trajectory upon emerging from the 
flow and therefore the polarization vectors started off mis¬ 
aligned with one another. Secondly, the parallel transport 
of the polarization vector in heavily curved space-time fur¬ 
ther changes the orientation. These effects are, as expected, 
stronger for photons which passed closer to the BH, including 
photons from the back of the flow that emerged from rela¬ 
tively large radii but needed to pass very close to the BH to 
reach the observer. We also see an asymmetry in this effect 
which is directly due to the frame dragging effect. The ori¬ 
entation of this vector oscillates as a function of precession 
angle because the orientation of the flow itself is oscillating. 
This is diluted by GR effects but not entirely. 

Figure [6] shows the (normalised) integrated flux, polariza¬ 
tion degree and polarization angle plotted against QPO phase 
for this run of the model, with solid and dashed lines rep¬ 
resenting respectively the full calculation and the simplified 
Newtonian approximation. We see that relativistic effects 
wash out variability in the flux and the polarization signature 
as well as reducing the average observed polarization degree 
and angle. The amplitude of variability is reduced mainly 
through light bending, which allows us to see the back of the 
flow even when the angle between the flow spin axis and the 
line of sight is large. The overall polarization degree is lower 
for the full calculation because the polarization vectors ob¬ 
served from different regions of the flow are not aligned, and 
therefore do not add together completely constructively as in 
the Newtonian approximation. GR influences the polarization 
angle because the polarization vectors of emergent photons 
appear to be bent around the BH. The vectors in Figure [5] to 
the left of the BH are forced clockwise and those to the right 
are forced anti-clockwise. This is a more pronounced effect 
on the left hand side of the BH because the ‘critical point’, 
where only light rays that are emitted perpendicular to the 
flow mid-plane can reach the observer, is situated here (due to 
frame dragging). At this point, the observer sees the polariza¬ 
tion vector rotate, mainly due to spe cial relativistic aberration 
(see e.g. Figure 3 in lDovciakll201C)l) . 

For the full calculation, the flux has a fractional rms of 
10.0% and 2.24% in the first and second harmonics respec¬ 
tively. This is representative of QPOs in high inclination 
sources (e.g. GRS 1915+105: lYan et al.llMS ) for the first 
harmonic but the amplitude of the second harmonic is slightly 
lower than is typically observed. For the polarization signa¬ 
ture, we instead consider the absolute rms, since this is rel¬ 
evant for detection. This is respectively 1.0% and 0.3% for 
the first and second harmonics of the polarization degree and 
respectively 3.0° and 0.2° for the first and second harmonics 
of the polarization angle. Also, the polarization degree lags 
the flux by 147.3° for the first harmonic but the polarization 
angle leads the flux by 125.4°. To estimate the importance 
of the secondary image, we additionally perform a calcula¬ 
tion whereby rays that pass under the disk plane and hit the 
underside of the flow are assumed to be blocked. For this al¬ 
ternative calculation, the flux has a slightly higher fractional 
rms of 13.3% and 3.75% for the first and second harmonics 
respectively. This is because the secondary image appears to 
have a roughly constant shape to the observer and so serves to 
wash out variability. The polarization properties are remark- 
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Figure 8. Images of the flow for three precession angles for the model with i = 30°, <I> = 180°, [3 = 10° and h/r = 0.1. The left hand side 
pictures blueshift and the right hand side pictures flux with the polarization vector overlaid, normalised to the maximum observed polarization degree. 
The three precession angles pictured, in units of cycles, are lj = 0, 0.3125 and 0.625 from top to bottom respectively. Full movies corresponding to 
these images can be found at http://figshare.com/articles/Polarization_modulation_gifs/1351920 . Individual movies can alternatively be found on YouTube: 
www.youtube.com/watch?v=TSe—iXofu8&feature=youtu.be (blueshift) and www.youtube.com/watch?v=GjlIRfkor_s&feature=youtu.be (flux and polariza¬ 
tion). 
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ably similar for the alternative calculation: the polarization 
degree has an rms of 1.0% and 0.3% for the first two harmon¬ 
ics and the angle has an rms of 3.1° and 0.1°. 

Figure[8]shows images of the low inclination model, for the 
same three precession angles as previously depicted in Fig¬ 
ure [5] For both the flux (right) and blueshift (left), we still 
see Doppler shifts due to rapid rotation, gravitational redshift 
and the effects of light bending, however all these effects are 
less pronounced since the component of velocity in the line 
of sight is lower and also less photons need to pass close to 
the BH in order to reach the observer. The full movies for 
all of the images shown in Figures [5] and [8] can be found at 
http://figshare.com/articles/Polarization modulation gifs/1351920 

(additionally, YouTube links for individual movies are given 
in the Figure captions). 

Figure [7] shows the flux, polarization degree and polariza¬ 
tion angle plotted against QPO phase for this model, again 
with the results of the simple Newtonian calculation repre¬ 
sented with dashed lines. Here, the flux is less affected by 
light bending, since the angle between the flow spin axis and 
the line of sight never gets particularly large. GR effects only 
reduce the mean and rms of the polarization degree and angle, 
with the phase of the oscillations only slightly modified from 
the Newtonian approximation. The amplitude of flux variabil¬ 
ity is lower here than in the high inclination case, although the 
inclination dependence of amplitude is far more pronounced 
for the Newtonian approximation than for the full calculation. 
This is again down to light bending. The absolute rms of the 
polarization degree is comparable to the high inclination case, 
but the fractional variation in polarization degree is actually 
significantly larger here than for the high inclination model. 
This can be understood from Figure [4] which shows that the 
emergent polarization degree is a far steeper function of 
for large p e . However, the mean polarization degree is far 
lower here than for the high inclination angle, and these two 
considerations happen to balance. The amplitude of the po¬ 
larization angle oscillation is greater for the low inclination 
model, which is purely a geometrical effect. It is clear from 
Figure [7] that the polarization degree and flux are nearly in 
phase for this example, in contrast to the results for the high 
inclination model. We find that p lags F by 20.3° for the fun¬ 
damental and x l a g s F by 113.8°. Since we see no secondary 
image from the underside of the flow for this viewing angle, 
assuming the material between the disk and flow to be opti¬ 
cally thick does not change our results. 

4.2. Parameter study 

We now consider the full range of viewing angles for the 
specific geometrical setup described at the start of this sec¬ 
tion. The corresponding parameter exploration in VPI13 con¬ 
sidered viewing angles ranging from 0° < i < 90° and 0° < 
$ < 180°. However, our use of the Kerr metric introduces 
a subtle asymmetry (which disappears for a = 0 of course) 
meaning that we must explore the range 0° < $ < 360° to 
be exhaustive. We use a reduced resolution of 80 x 80 pixels 
and consider 16 precession angles. When we test for the spe¬ 
cific high and low inclination models considered above, we 
find that this reduced resolution provides a very good approx¬ 
imation to the high resolution run. 

In Figure [9j we use colors and contours to plot four dif¬ 
ferent quantities for the entire range of viewing angles. Each 
plot has it’s own logarithmic scale. The top-left plot shows the 
fractional rms amplitude of the flux modulation, op / (F), for 
the first harmonic only. The contours in this plot are labeled as 


percentages and labels always face in the uphill direction. We 
see that the amplitude of the QPO increases with i as we may 
expect from t he previous subsection. This is consistent with 
obser vations (ISchnittman et alJ2006allMottall2014l : lHeil et afl 
120151) . We note that this rise of amplitude with inclination 
angle would be even more pronounced if we had assumed 
secondary images to be blocked by optically thick material, 
since these wash out variability for high viewer inclinations. 
The plot (and all subsequent plots) is nearly symmetric about 
$ = 180° but a subtle asymmetry is introduced by the as¬ 
sumption of high spin. The top-right plot shows the absolute 
rms amplitude of the polarization degree modulation a p , again 
for the fundamental only, with contours labeled as percent¬ 
ages. This peaks at i ~ 60 — 70°. The bottom-right plot shows 
mean polarization degree (p), again with contours labeled as 
percentages. We see that this increases with i. For small i, the 
mean polarization degree is very low indeed, which leads to 
the fractional rms amplitude of the polarization degree modu¬ 
lation becoming very large for low inclinations. However, it is 
likely the absolute rms amplitude that is relevant to detection. 

The bottom-left plot shows the absolute rms variability of 
polarization angle o x , with contours now labelled in degrees, 
again for the first harmonic only. Note that the polarization 
angle is only defined on an interval of 180°, since upward 
polarization is indistinguishable from downward polarization. 
This can produce phase wrapping: e.g. if x oscillates between 
a minimum and maximum of 80° and 100° but is defined on 
the interval —90° to 90°, the measured x will jump from 90° 
to —90° at some point during the QPO cycle, producing a 
spuriously large rms measurement. To avoid this, for every 
combination of i and <1>, we define x for the first QPO phase 
on the interval —90° to 90°. For each subsequent value of 
QPO phase, we first calculate x °n the same interval, but also 
try adding 0, 180° and —180°, and choose the interval that 
minimises the difference to the previous measurement of X- 
We see that the amplitude decreases with inclination angle, 
but is above 1° for most of parameter space. 

In Figure [TO] (left), we plot the phase lag between polar¬ 
ization degree and flux against i and <f>. Here, positive lags 
mean that p lags F and we again consider only the fundamen¬ 
tal. We see that p lags F for nearly all viewing angles. In the 
right panel, we instead plot the phase lag between x and F. 
Here we see phase wrapping at i ~ 60° which occurs because 
the lag between the two functions is only defined on the inter¬ 
val — 180° to 180°. The magnitude of the lag is large for most 
of parameter space. This can be understood by looking at the 
images in Figures (0 and ©. The peak in x for our coordi¬ 
nate system occurs approximately when we see the flow spin 
axis tilted the furthest to the left. The flux rarely peaks close 
to this phase in the precession cycle. If a high enough sig¬ 
nal to noise can be achieved to observe these lags, they may 
provide a powerful extra diagnostic. 

5. DISCUSSION 

We have calculated the polarization signature predicted by 
the precessing inner flow model for low frequency QPOs. We 
find that the polarization degree and angle are expected to be 
modulated on the QPO period. 

5.1. Assumptions 

In our analysis, we calculate the GR effects very accu¬ 
rately but make a number of simplifying assumptions about 
the properties of the inner accretion flow. We use analyti¬ 
cal parameterizations for the angular dependence of intensity 
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Flux: fractional rms (%) 


Polarization degree: absolute rms (%) 




i (degrees) i (degrees) 

Figure 9. Top-Left: Fractional rms in the first harmonic of the flux modulation plotted as a function of the two viewing angles, i and <f>, with contours labelled 
as percentages. Top-Right: Absolute rms in the first harmonic of the polarization degree modulation, with contours again labelled as percentages. Bottom-Right: 
Mean polarization degree, with contours once again labelled as percentages. Bottom-Left: Absolute rms in the first harmonic of the polarization angle modulation, 
with contours now labelled in degrees. The colors in each plot follow separate logarithmic scales. We see that the amplitude of the flux and polarization degree 
modulations increase with inclination angle, i, whereas the amplitude of the polarization angle modulation reduces with i. Note the slight asymmetry in all plots 
around = 180°, which occurs because a/0. 
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Polarization degree lagging flux (degrees) 


Polarization angle lagging flux (degrees) 




0 20 40 60 80 
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Figure 10. Left: Phase lag in degrees between the polarization degree and the flux, with positive lag meaning that p lags F. Here, we only plot for the 
fundamental and color scales linearly from —180° (darkest) to 180° (lightest). We see that the model predicts the oscillation in polarization degree to lag the 
first harmonic of flux for almost the full range of viewing angles. Right: Phase lag between polarization angle and flux. The model predicts the magnitude of the 
phase lag between x and F to be large for most of parameter space. 


Polarization degree: required count rate (c/s) 


Polarization angle: required count rate (c/s) 
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Figure 11. Required count rate in order to detect a modulation in polarization degree (left) and angle (right). We assume a 200 ks exposure and a modulation 
factor of // = 0.5. See text for more details. 
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and polarization degree of radiation emerging from the flow. 
Since these parameterizations are based on the calculations 
of ST85, this is a reasonable assumption but there is scope 
to extend this work in future. Our parameterization is only 
valid for one value of optical depth (r = 1) and so we cannot 
explore the dependence of our results on this parameter. In 
addition, we assume that the angular dependence of the flow 
emissivity can be separated from the radial and energy depen¬ 
dencies (Equation [9]). In reality, it is likely that the shape of 
the emitted spectrum depends on both radius and viewing an¬ 
gle. Since the outer part of the flow is illuminated by a greater 
flux of cool disk photons, the spectrum should be softer for 
larger r, g iving rise_to the observed time lags dKotov et al.l 
1200 It llngram & van der Klisl 120131) . There is now observa¬ 
tional evidence that th e spec tral shape of the flow depends on 
viewing angle dHeil et al.12015ilTngram & van der Klisi2015l) . 
Similarly, the angular dependence of polarization degree will 
likely depend on energy of the emitted photon, and also per¬ 
haps the radius it was emitted from. We also expect the flow 
optical depth, and therefore mean polarization degree, to de¬ 
pend on truncation radius / spectral state. For a low luminos¬ 
ity hard state we expect r to drop below unity and polarization 
degree to reach ~ 50% (see lViironen & Poutanenll2004l) . Re¬ 
placing our simple parameterization with a full Monte Carlo 
simulation in future will be a significant improvement upon 
this work. 

We also effectively assume that the seed photon luminos¬ 
ity stays constant during a QPO cycle. This is a good as¬ 
sumption if the seed photon luminosit y is com pletely domi¬ 
nated by internally generated photons (IVeledina et al.ll201 1 ak 
VPI13). However, the luminosity of disk photons incident 
on the flow will vary as the misalignment between the disk 
and flow changes over the course of a precession cycle. A 
calculation of disk seed photon luminosity as a function of 
precession angle has yet to be performed taking GR effects 
fully into account. The relative importance of internally gen¬ 
erated seed photons will increase with truncation radius. It 
should be possible to put constraints on such an evolution of 
the seed photon origin with X-ray polarimetry, since polariza¬ 
tion degree depends strongly on scatte ring order for disk seed 
photons (see Fig. 2 in lViironen & Poutanenll2004l) . 

5.2. Detection 

It may be possible to detect the predicted polarization de¬ 
gree and angle modulation in the near future with a dedicated 
polarization satellite mission. Here we present a simple cal¬ 
culation to roughly assess detectability. X-ray polarimeters 
generally use either Thomson scattering (e.g. the Polarization 
Spectroscopic Telescope Array, PolSTAR) or the photoelectric 
effect (e.g. Gravity & Extreme Magnetism SMEX\ GEMS). 
Scattering polarimeters measure the landing position of scat¬ 
tered photons, which travel preferentially in the direction per¬ 
pendicular to their electric field vector (e.g. IGuo et al.ll2013h . 
Photoelectric effect polarimeters track the direction of pho¬ 
toelectrons, which are preferentially emitted in the direction 
parallel to the incoming electric field vector (e.g. iBlack et aD 
120101) . In either case, an estimate for the polarization angle 
of each incident photon, Xk (he- for the k th photon), can 
be recorded. When many photons are incident on the detec¬ 
tor, a polarized si gnal will exhibit a sinus o idal modulation 
with distribution ( Krawczvnski et all 1201 U iFei et al.l 1 1 997b 
iKislat et alMll) : 

f{Xk) oc l+p/rcos[2(xfc-x)], (30) 


where \ i s the ‘true’ polarization angle of the source, p is 
the ‘true’ polarization degree and p is the modulation factor. 
The performance of a polarimeter can be characterised by the 
modulation factor, since it governs the distribution of a 100% 
polarized signal. When C photon s are i ncident on the detec¬ 
tor, the measurement error on p is (IKislat et al.ll2014l) 


dp 


2 /p 1 — p 2 

C-l 


(31) 


when C is sufficiently large to be in a Gaussian regime (note, 
p is represented here as a fraction rather than a percentage). 
The measurement error on \ (in radians) is 


dx 


1 

ppsj2(C - 1) 


(32) 


Detection of a ~ Is QPO requires a time resolution < 0.1s. 
Very few photons can likely be collected in such a short time 
and therefore a high time resolution time series of individual 
Stokes parameters would be very noisy indeed. It is , however, 
possib le to stack into QPO phase bins dTomsick & Kaaretl 
[20011) . resulting in one folded time series. If we stack into 
TV phase bins, the number of counts detected per phase bin 
is simply C = RT /TV, where R is the mean count rate and 
T is the total exposure time. We can estimate that detection 
of a QPO in the polarization degree requires a measurement 
error in each phase bin around one fifth of the absolute rms 
variability amplitude, dp ~ o v j 5. Thus, the mean count rate 


we must be able to detect for the observation is 


R ■ 


25 TV 


T<J l IP 




(33) 


We can estimate the count rate required to detect a modulation 
in polarization angle in a similar manner. In this case 


R ~ 4.1 x 10" 


TV 


T(p) 2 p 2 o 2 ’ 


(34) 


where o x is in degrees. In Figure[TTl we show these required 
count rates for all of parameter space assuming TV = 8 phase 
bins, an exposure of T = 200 ks and a modulation factor of 
p = 0.5. The plots on the left and right are for polarization 
degree and angle respectively. We see that to detect either 
modulations at all, a polarimeter must have the sensitivity to 
measure a count rate above 40 c/s, and a count rate of ~ 60 c/s 
opens up a reasonably large fraction of parameter space (i ~ 
40° — 75°) in which both modulations can be detected. The 
same folding method can be used to measure lags between the 
polarization properties and the flux, although this may require 
a more challenging sensitivity. 

A number of missions proposed for NASA’s Small Explorer 
Program (SMEX) should be able to detect count rates > 40 
c/s. PolSTAR (the satellite incarnation of the balloon experi¬ 
ment X-Ca/i7>«r; |Guoetai]20T3l) is sensitive to the 10 — 20 
keV energy range in which the Comptonized spectrum dom¬ 
inates, and would therefore be ideal for this application. The 
soft X-ray missions GEMS (IBlack et al.ll2010l) and the Imag¬ 
ing X-ray Polarimetry Explorer (/XPEf lWeisskopf et al.l2008l) 
should also be able to measure a similar count rate, how¬ 
ever contribution to the flux from the constant disk compo¬ 
nent at softer energies will dilute the variability amplitude of 
the polarization signature, making detection more challenging 
(rms ss rmsfl 0 w[l — £disk], where Xdisk is the fraction of the 
flux contributed by the disk). Disk dilution is not a problem 
in the hard state, but here the flux is lower and the QPOs are 
generally not as coherent as in the HIMS. 
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5.3. Implications 

Detection of a polarization modulation on the QPO fre¬ 
quency would have strong implications. First of all, this 
would confirm that the QPO is indeed a geometric effect, 
as is strongly hinted in the literature at the momen t dMottal 
l2014t iHeil et all 1201 5t llngram & van der Klisll2015l) . It will 
also provide a strong test for the Lense-Thirring QPO model. 
This is fairly profound in itself, since Lense-Thirring preces¬ 
sion has never been unambiguously observed in the strong 
field regime, but this also has implications with regard to mak¬ 
ing spin measurements_jCurrent spin meas urements from disk 
spectr al fitting (e.g. lKolehmainen & Donel201CtfSteiner et al.l 
1201 ll) assume that the BH and binary spin axes are aligned, 
in contradiction of the Lense-Thirring QPO model. Care¬ 
ful modelling of QPO properties can be used to estimate the 
misalignment angle /3 to improve the spectroscopic measure¬ 
ments and also the very measurement of the QPO frequency 
itself gives an orthogo nal spin estimate if high frequency 
QPOs are also present dMotta et al .1120141 : llngram & Mottal 

BoTl - 

6. CONCLUSIONS 

We find that the polarization signature emitted from a trun¬ 
cated disk / precessing inner flow geometry oscillates on the 
QPO frequency. The modulation in polarization degree has 
an absolute rms which increases gradually with viewer incli¬ 
nation angle, peaking at 1.5% for i ~ 60°. In contrast, the 
absolute rms of the polarization angle modulation is higher 
for lower inclination angles. For polarization degree, this in¬ 
clination dependence is mainly due to the assumed angular 
dependence of polarization degree and emissivity for Comp¬ 
ton scattering, which we parameterize to agree with the calcu¬ 
lations of ST85. Although our parameterization is only appro¬ 
priate for an optical depth of r = 1, we note that the results of 
ST85 are qualitatively similar for optical depths up to r ~ 2.5. 
The inclination dependence of the polarization angle depends 
mainly on the precessing geometry and GR effects, and so is 
robust to our assumptions about the Compton scattering pro¬ 
cess. Our calculations here only consider one specific flow 
geometry, and we explore the full range of viewing angles. In 
future, we will also explore the effects of changing parameters 
such as the truncation radius and spin. In particular, we note 
that assuming a larger misalignment angle, /3, will increase 
the predicted amplitudes of all modulations. 

We find through a rough calculation that, in order to detect 
this effect for a reasonable fraction of parameter space, an 
X-ray polarimeter will need to detect a 10 — 20 keV count 
rate of > 60 c/s from a bright object displaying QPOs. The 
current generation of proposed X-ray polarimetery missions 
will likely fill this requirement. In particular, NASA’s PolSTAR 
is suited to this application due to its sensitivity to hard X- 
rays. 
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the Academy of Finland for support (grant 268740). We thank 
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APPENDIX 


A: THE KERR SPACE-TIME 

The non-zero entries of the Kerr metric can be expressed in Boyer-Lindquist coordinates as 

—2 ar sin 2 8 E 

) = JT ) 9<t>t = 9tcf> ) 9rr = > 988 = 



A sin 2 8 

E 


(Al) 


where E = r 2 + a 2 cos 2 8, A = r 2 — 2r + a 2 and A = (V 2 + a 2 ) 2 — Aa 2 sin 2 8. Th e 4-momentum of photons travelling in the 


Kerr space-time can be expressed as (!Carterlll968l : lMisner et al.fl 9731 : lDovciakll20()4l) 

p t = E _1 [ a(l — asin 2 8) + (r 2 + a 2 )(r 2 + a 2 — al)/ A] (A2) 

p r = i? sgn E —1 [(r 2 + a 2 - al) 2 - A[(7 - a) 2 + q 2 ]] ^ (A3) 

p e = —0 S g n E _1 [q 2 — cot 2 0(l 2 — a 2 sin 2 0)] 1 ^ 2 (A4) 

p^ = E _1 [(/sin 2 8 — a + a{r 2 + a 2 — a/)/A] , (A5) 


where Carter’s constants of motion, l and q 2 are given in the main text and the sign of the radial and polar coordinates is denoted 
by f? sgn and 0 sgn . Both of these are positive for the case of an even number of turning points between the viewer and the emission 
point and negative for an odd number. 


B: THE FLOW NORMAL 


According to the equivalence principle, we can always define a ‘free falling laboratory’ frame in which GR reduces to Special 
Relativity. Mathematically, this can be achieved by defining a tetrad of orthonormal unit 4-vectors. We require the flow normal riJ' 
to be part of such an orthonormal tetrad in order to calculate the emission angle from it. We calculate the r, 8 and </> components 
of nJ 1 by transforming Zf into Boyer-Lindq uist coo rdi nates ( see Appendix 0. The time-like unit vector in the tetrad is simply 
u ,J ", i.e. the instantaneous rest frame (e.g. lKrawczvnskill20l 2t IWil kins & F abian 2012). The time-like component of n 1 ' can thus 
be calculated by setting = 0 to give 


n t = Zf9disU L ' 

9tfiU p 


(Bl) 


where d = r, 8, </> but Greek letters run from t to 0 as usual. Finally, we normalise to ensure iiJ'rif, = 1. 


C: COORDINATE TRANSFORMS 

We can convert from the BH angles 8 and (j) to the flow angles 8f and Of using the formulae 



cos 8f = ? • Zf 

(Cl) 


, r • yf 

tan <j)f =-. 

r • Xf 

(C2) 

cos Of = 

sin 8 sin /3 cos(w — ujq — 4>) + cos 8 cos /3 

(C3) 

tan 4>f = 

sin 8 sin(o; — loq — (f>) 

(C4) 

cos 8 sin /? — sin 8 cos /3 cos(u; — ujq — (/>)' 


This gives 
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The reciprocal conversion is 


and 


tan <b = 


cos 9 = cos Of cos/3 + sin#/ cos 0/ sin/3 

— sin0/ [cos /? sin© — wo) cos 0/ + cos(w — wo) sin 0/] + cos Of sin /3 sin(w — wo) 
sin© [— cos/3 cos(w — wq) cos 0/ + sin(w — wq) sin 0/] + cos 0/ sin/3 cos© — wq) 


(C5) 

(C6) 


Additionally, since Zb = Zf (w = 7r), we can convert to and from binary angles using the above formulae by setting w = tt (i.e. 
replacing subscript f with subscript b and replacing w with 7r). 

Our calculation also requires the conversion of Cartesian coordinates to Boyer-Lindquist coordinates. The two are related as 


x = Vr 2 + a 2 sin 0 cos 0 
y = Vr 2 + a 2 sin 0 sin 0 
z = rcos9. 


(C7) 


To convert the coordinates of a 3-vector, A, from a Cartesian to a Boyer-Lindquist representation, we simply apply the general 
formula 

dx b ' 


A b ' = 


dx a 


■A a 


where b' can take the values 


1 and 0 and a can take the values x, y and z. Combining Equations (lC7b gives 

o l r 


where p 2 = x 2 
expressed as 


2 L‘ 


p 2 - a 2 + \J(p 2 — a 2 ) 2 + 4a 2 z 2 


(C8) 


(C9) 


z 2 , along with the familiar expressions c:os 0 = zjr and tan 0 = y/x. The differentials of r can be 


and for i 


dr 

dx 

dr 

dy 

dr 

dz 

d6 

dx 

d0_ 

dy 

d6 

dz 


= \ 1 + 


a N 2 sin 6 cos <fi 
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1 + 
cos 8 


a\ 2 sin 9 sin (f> 
~2 

■a 2 (H 


© 


1 + 


1 + 


— a 2 cos 2 0" 1 


— a 2 cos 2 9 


1 + 


■ sin 2 9) 


(CIO) 


= \ 1 + 


© 


1 + 
r sin 9 


© 


2 cos 9 cos </> 
2r 

! cos 9 sin 4> 


1 + 


r 2 — a 2 cos 2 


2 r 


1 + 


(Cll) 


Note that these expressions, as expected, reduce to the relations for spherical polar coordinates in the Schwarzschild limit. The 
case of sin 9 = 0 must be treated separately, although this is rather stra ight forward. The differentials of 0 are the same as the case 
of spherical polars. These can then be substituted into Equation ( 1C81 in order to conve rt the vector. The reciprocal conversion, 
from Boyer-Lindquist to Cartesian, is far simpler requiring the differentials of Equations (lC7b with respect to the Boyer-Lindquist 
coordinates. 


































